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Abstract 

A one phase Stefan problem in nonlinear conduction is considered. The problem is 
shown to admit a unique solution for small times. An exact solution is obtained which 
is a travelling front moving with constant speed. 



One and two phase Stefan problems for the linear heat equation have been the subject 



1 Introduction 

> 
00 

of many studies in the past [1, 2]. Indeed these problems have a great physical relevance 
since they provide a mathematical model for the processes of phase changes [3, 4]. 

The boundary between the two phases is a free boundary: its motion has to be deter- 
mined as part of the solution. 

More recently the previous analysis was extended to nonlinear diffusion models. In [5, 6] 
Q_i' exact solutions were found in parametric form for a class of Stefan problems in nonlinear 

heat conduction. 

Moreover one and two-phase Stefan problems for the Burgers equation were solved 
in [7, 8] and explicit travelling wave solutions were obtained. 

It is the aim of this paper to formulate and solve a one-phase Stefan problem for the 
nonlinear heat equation: 



i) 



= *om), * >0 > ( L1 ) 

on the semiinfinite domain x G (— oo, s(t)) characterized by the following set of initial and 
boundary data 

$(x,0) =tio(x), xe(-oo,b) (1.2a) 

and , i?oO) = Ih < 0, b > 0, 

0(-oo,t) =/3i >0, t>0, /3i>|0o(x)| > |/H 0x(-oo,t)=O, (1.2b) 

0(a(t),t) = fa, *>0, with s(0) = 6, (1.2c) 

Ms(t),t) = -s(t)- (l-2d) 
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Equation (1.1) a well know mathematical model for heat conduction in solid crystalline 
hydrogen for situations with one-dimensional spatial symmetry [9] . The one-phase Stefan 
problem (1.1)— (1.2) is associated to a phase change (fusion) in such a system. 

In the above relations (3±, 02 and b are constants; the unknown function s(t) describes 
the motion of the free boundary and has to be determined together with $(x, t). Moreover, 
equation (1.2d) is a condition on the flux at the free boundary, arising from heat balance 
(energy) considerations. 

In the following we assume $o(x) to be a continuously differentiable function of its 
argument. 

Our analysis is based on the approach followed in [7] for the solution of a one-phase 
Stefan problem for the Burgers equation. 

In the next Section we reduce the problem (1.1)— (1.2) to a nonlinear integral equation 
for the independent variable t. 

In Section 3 we prove existence and uniqueness of the solution for small intervals of 
time; in the last Section we show that the system admits an exact solution which travels 
with a constant velocity proportional to the velocity of the free boundary. 

Some details of the proof presented in the third Section are given in the Appendix. 

2 Linearization 

In order to linearize equation (1.1) we introduce the transformation 

ip(z,t) =ti(x,t), z = z(x,t), z x = -, z t = - l & x , (2.1) 

v 

whose compatibility z x t = zt x is easily proved via (1.1). Under this transformation equa- 
tion (1.1) is mapped into 

It = 1>zz (2.2) 

on the domain — oo < z < z(t) with z(t) = z(s(t),t) and z(0) = b. 

(2.2) is the linear heat equation for the dependent variable tp(z, t) with initial datum 
given by 

1>(z,0) = M*o) =Mx), (2.3a) 

where 

•"M (2 ' 3b> 

The boundary conditions (1.2b) and (1.2c) take now the form 

^(-oo,*) =0(-oo,*) =/?i, 

^ 2 (-oo, t) = 0a.(-oo, *)0(-oo, = (2.3c) 

and 



i>(z{t),t) = 02, 

*p z (z(t),t) = -fos(t). 



(2.3d) 
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The Stefan problem for the nonlinear equation (1.1) has then been mapped into a clas- 
sical Stefan problem for the heat equation (2.2) with initial datum (2.3a), characterized 
by the boundary conditions (2.3d) at the free boundary. 

We say ip{z,t), z(t) form a solution of the above Stefan problem for t < a, < a < oo, 
when: (i) tp(z,t) is a solution of (2.2) satisfying (2.3), it exists and is continuous together 
with its derivatives for — oo < z < z(t), < t < a; (ii) s(t) is a continuously differentiable 
function for < t < a. 

In the following we outline a method to prove the existence and uniqueness of the 
solution for small times, t < a. 

We first observe that by integrating the second relation in (2.3d) we get 



P2 Jo 

which in turn implies 

z(t) = h(t)-±- f v*(*(tVX 

P2 Jo 



with 



h(t) 



s(t) 1 



dx'. 



(2.4a) 



(2.4b) 



(2.4c) 



Next we turn our attention to the solution of (2.2). We introduce the fundamental 
kernel of the heat equation 

(*-0 2 " 



K(z-^t-r) 



1 



1 



20F V / F^7 



exp 



4(t-r) 



and integrate Green's identity for the heat equation 



d_( K d± ,9K\ d_ 



{Kil;) = 0, 



(2.5) 



(2.6) 



over the domain — oo < £ < z(t), e < r < r — e and let e — > 0. Using ij}{z{r), r) = and 
K(z — £, 0) = S(z — £)> we obtain 

i>(z,t)= f K(z-Z,t)MZ)d{i--^ f K(z-z(T),t-r)MHr),r)dT 

J-oo P2 Jo 



f K^z{t), 
Jo 



t — r)dr, 



(2.7) 



with z(t) and h(t) given by (2.4b) and (2.4c) respectively. 

In the right hand side of (2.7) ip z (z(t),t) is unknown; it is convenient to take the 
z-derivative of both sides in (2.7) and evaluate it as z — ► z(t)~ . 

By putting v{t) = ip z (z(t),t), we obtain: 



u(t) 



1 

W2 

■t 



-Mb)K(z(t) -b,t)+ f K(z{t) - e, i)V>o(£R 

J — 00 

-j- [* K z (z(t)-z(T),t-T)dT-fo f K T (z(T),t-T)dr] , (2.8a) 

P2 Jo Jo 
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with 

i r* 

z(t) = h(t) -— v(t)<It. (2.8b) 

P2 JO 

Thus the solution of the Stefan problem (2.2), (2.3a), (2.3d) has been reduced to the 
solution of the nonlinear integral equation (2.8a) and (2.8b) for the independent variable t. 

Once the existence and uniqueness of the function v{t) is established for < t < a, 
there follows via (2.7) the existence and uniqueness of ip(z, t) (and of i?(x, t)) for < t < a. 



3 Contraction Mapping 

In order to analyze existence properties of v(t) for < t < a, we denote by Sm the closed 
sphere \\v\\ < M in the Banach space of functions v(t) continuous for < t < a with the 
uniform norm ||z/|| = l.u.b. |^(t)|. On the sphere Sm define the transformation 



w = Tv, 



(3.1) 



where Tv coincides with the right hand side of (2.8a). We first prove that T is a mapping 
of Sm into itself. From (2.8b) we obtain 



W)\<Wt)\ + ^- 

\P2\ 



u{t)(1t 



< 



M 

W\ 



(3.2) 



where (2.4c) have been used. From (2.8b) we also get 



\m-z(t')\<\h(t)-h(t')\ + — 

\P2\ 



v(t)<It 



< 



M 

W\ 



1 



l+ M )lt-fl^llt-fl. 



We now turn our attention to the right hand side of (3.1). We first note that 
\K(z(t) -b,t)\ < B 2 a, 



(3.3) 



(3.4) 



where B 2 is an appropriate constant depending on b. Next we consider the integral terms 
in the right hand side of (3.1). We can write: 



K{ z {t) - i,t)^)di 



< W 



r~ b i i 



exp 



m - 

At 



< 



and 



/? 2 f K T {z{t),t-T)dT 

Jo 



1 1 



exp [-z 2 (t)/4t] < \P 2 \B 2 <J. 



(3.5a) 



(3.5b) 
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Moreover we get 



-j- f K z {z{t)-z{t'),t-t')v{t')dt' 

P2 Jo 

M 1 /•* \z(t) -z(t')\ , / M 1 



2|/3 2 |v^yo — f)i 



W < 



(3.5c) 



where (3.3) has been used. 

We now use the inequality ( 1 



2|A| 



> | and define M as M = 2Ai + 1. When (3.1) 



is used together with (3.4) and (3.5) we get 
IMI < M 



(3.6) 



provided we choose cr = min(cri, (T2) with c7i : (I/^I + IV'oWD-^o'i < j and 02 '■ MB\^fo2 < 
^jy^ ; thus the mapping is closed. 

Next we wish to prove that T is a contraction; i.e. given two solutions of (3.1) with 
\\v — v\\ = 5, it follows that \\tu — ti>\\ = $8 with < $ < 1. Using (2.8b) we have for small 
enough 5 



1 



\z(t) ~ Ht')\ < T^St 



s(t) 1 



dx' 



< 



1 

W\ 



St 1 + 



1 

W\ 



= B 3 St, 



(3.7) 



where (2.4c) has also been used. 

Similar estimates hold for z(t), which will be useful in the following. From (2.8b) we 
see that z(t) is bounded 



1 



M 



Kt)\ < \Ht)\ + T7nW(t)\ < 7^7 1 + T5-T = 5 3M, 



1 



moreover it is 



^From (3.1) we now write 



w — w 



-1 4 



1 + 



(3.8a) 



(3.8b) 



(3.9a) 



with 



(3.9b) 
(3.9c) 



H 1= Mb) [K{z(t) -b,t)- K(z(t) - b, t)] , 

h 2 = f $,(0 [Km) -z,t)- K(m - e, t)] 

H 3 = ±- f [K z (z(t] - 1(0, * - 0K0 - ^*(*(f) - *(0>* - 0*(0] < ( 3 - 9d ) 

P2 Jo 

H A = {3 2 f [K t ,{z(t),t-t')-K tl {z(t),t-t')}dt'. (3.9e) 

JO 
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First we estimate H±. We use the mean value theorem together with (3.7) and (3.2) in 
the right hand side of (3.9b); we get 



|#i| < \\M\-^T=B 1 B z 5t<^ i B 1 B^5 = B^5. 
The estimate of H 2 in (3.9c) is obtained by writing 



(3.10) 



m<W[ 

/7T 



dy 



(3.11a) 



with y = ^yf- and y = z ^^- fe ; we then obtain 

\H 2 \ < Wl^im-Kt)] <^B 3 ^S = B 5 ^5, (3.11b) 

where (3.7) has been used. 

The estimate of H% is somewhat more cumbersome; a detailed analysis is given in the 
Appendix (see (A.1)-(A.6)). 

There obtains 

\H Z \ < B 6 ^5 (3.12) 

where Bq is defined in (A. 7). 

We finally turn our attention to the estimate of H4 in (3.9e). When (2.5) is used, we 
get from the integral in the right hand side of (3.9e) 

\H 4 \ < -^L |exp(-l 2 (t)/4t) -exp(-I 2 (t)/4t)|< -j^L B^z^/ad = B 7 y/aS, (3.13) 

2\/TTt 4V7T 

where use of the mean value theorem together with (3.2) and (3.7) has been made. 
/,From (3.9a) we now write 

l--H< TTW gl^l, 

which in turn implies, when we combine together the estimates (3.10)-(3.13): 

7 



\w — w 

6 , 1 



< V^J2 B i = V^^85 (3.14) 



thus we conclude that if a satisfies a < min(<7i, 02, 03) where 

V^B$ < a (3.15) 

it follows that T is a contraction operator on Sm, which admits a unique fixed point 
v = Tv in Sm for < t < a. 

We have then proven the existence and uniqueness of the solution of the integral equa- 
tion (2.8a) for a small interval of time. 
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4 A Particular Solution 

We now turn our attention to a particular solution of the Stefan problem (1.1), (1.2). 
Namely, we consider a moving front solution of equation (2.2) 

#M)=/?i(l-e- V( *- Vt) ), (4.1a) 

with 

V < 0, (4.1b) 

which is travelling to the left with constant speed V and is compatible with the boundary 
conditions (2.3c). We now impose on (4.1a) the Stefan boundary conditions (2.3d): the 
first one implies 

z = b + Vt, (4.2a) 
which in turn gives 

j = F =-M i+ f)- <42b> 

The boundary z(t) and the front solution (4.1a) are then both moving to the left with 
the same constant velocity. 

When we next use the second boundary condition (2.3d), keeping into account the first 
one, we obtain 

s{t)=®^M Vt (4.3) 

which shows that the moving boundary s(t) of the Stefan problem (1.1)— (1.2) is moving 
to the left with constant speed s = aV, a = ^j^- > 1. 

Finally, the solution of the one-phase Stefan problem for the nonlinear heat equa- 
tion (1.1) is given by 

( 4 - 4 ) 



dx r 

where, in virtue of (2.1), z(x,t) solves 



x 



= r^z',t)dz', (4.5) 
Jo 



with ip(z,i) given by (4.1a) and the speed V specified by (4.2b). 

We emphasize that the above solution is a very special solution of the Stefan prob- 
lem (1.1), (1.2). Indeed, it corresponds to particular case when the nonlinear integral 
equation (2.8) reduces to a linear integral equation of Volterra type in t, as implied by 
substituting back (4.2a) into (2.8a). 
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Appendix 



In order to estimate H 3 , starting from (3.9d) we write 
1 



H3 = T 

P2 



1 efrp(o * ( *|_yV (f(t) - i(o, t - 

^.u,At)-z(t') 



+ J dt'v(t') ^> t _J > K(z(t)-z(t'),t-t') 



with K(z — £, t — t) given by (2.5). 
Next, we put 

H 3 = ±-(J 1 + J2 + J3), 

P2 



Ji = - jT { W) - K0) ^_f° *W) - S(0, f - 0} , 



J 2 



j l dt'iu{t') 



f (t) - z(t') z(t) - z(t') 



t-t' 



t-t' 



K(i(t)-i(t'),t-t')\, 



■h - - f dt'[u{t') ~ zit \ f' ] K(z(t) - z(t'),t- t>) 







1 — exp 



t-t' 



By using (3.3) we estimate J\\ 



4(t-tf) 



l(t)-"(0 



t-t' 



v^-t' Vv 7 ^ 



The estimate of J2 uses (3.8b) and the mean value theorem: 



f (t) - f (f) f (t) - z(f) 



t-t' 



Finally, for the estimate of J3, we call 

_ (i(t) - mf - m - -z(t')? 

W " 4(t - t') 

_ p(t) - 1(f)) - (z(t) - z(t'))][(z(t) - z(t')) - (f (t) - l(t'))] 

4(t-t') 

Then from (3.3) and (3.7) we have 
\Q\< I ^ rV ^B 1 B 3 \t-t'\a5, 



(A.l) 



(A.2a) 
(A.2b) 
(A.2c) 



(A.2d) 



(A.3) 



(A.4) 



(A.5a) 
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which gives 

\Q\ < B^aS. (A.5b) 
We also note that we can estimate Q via (3.3): 

\Q\ < W) ~ ^')) 2 + (*(*) - *(0) 2 ] 

< X -B\\t - t'\ < \b\o < Bf, | < 1. (A.5c) 
Using |l — e _< 3| < |Q|e^l and (3.3) we estimate J3: 



m - m 



[ \l-e-Q\dt' 



< M= BfB^a^S < ^M Bl B 3 e B 'V^j 5, (A.6) 

where the definition of a following (3.6) has also been used. 

Combining the estimates of Ji, J2 and J3, we have from (A2.a) 

\P2\ Vv 71 " v 71 " \P2\ 4 y 
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